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I. Introduction 



A. The Weyl-Lanczos equations and the Lanczos Tensor Wave 
Equation in 4 Dimensions 

Lanczos ^2] generated the spacetime Weyl conformal tensor Cabcd from a 
tensor potential Labc by covariant differentiation such that Cabcd is given by 

Cabcd = Labc;d — Labd;c + Lcda;b — Lcdb;a + QbcL^ad) + 9adL(^hc) ~ 9bdL(^ac) 
—gacL{hd) + -^L'^'^ rn-^s{9ac9bd — 9ad9bc) , (1) 

where = La^ d;s ~ I^a'^s;d ^^id denotes covariant differentiation. We 
call the Weyl-Lanczos equations. The index symmetries of the Lanczos 
tensor Labc have to match the symmetries of and so it is usual to impose 

Labc = L[ab]c ) (2) 

and 

L[abc] = (3) 
and the trace free (gauge) condition 

La% = . (4) 

The spacetime Weyl Lanczos equations (Q) can also be expressed in the more 
compact form 

Cabcd = L[ab][c;d] + L[cd][a;b] "* L*[ab][c;d] ~* L*cd][a;b] i^) 

as done in [TH], where denotes covariant differentiation. The algebraic 
equations (jSJ, (jU leave us with only 16 independent components for the 
Labc- If we then introduce the differential gauge conditions 

Lab% = , (6) 
we can simplify (0) considerably to get 

Cabcd = Labc;d ~ L abd\c'^ L cda\b ~ Lcdb;a~ 9bcLad-s ~ 9adLlc-s ~^9bdLacg~^9acLld-s ■ 

(7) 

Theoretically, we could completely solve the 6 differential gauge conditions 
for 6 further components and have 10 components for 10 independent space- 
time Weyl-Lanczos equations. But this approach does not exhibit the most 
general solution. We note that equations © , Q constitute a system of lin- 
ear first-order partial differential equations in 4 dimensions which can easily 
be rewritten as an exterior differential system EDS in involution. This the- 
ory only applies in 4 dimensions as the Weyl Lanczos problem in 2 and 3 
dimensions does not exist and for 5 dimensions, we expect extra conditions 
to apply as we shall point out in a later paper. It was already shown in 
^ that the above 4-dimensional Weyl-Lanczos equations consist of an ex- 
terior differential system (EDS) in involution with respect to the spacetime 
variables. One can obtain the same results by applying the corresponding 
Janet-Riquier theory |H]. 

From the Weyl-Lanczos problem it is possible to generate a tensor wave 
equation for the (spacetime) Lanczos potential from which the Penrose wave 
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equation for the Weyl tensor Cabcd can be derived [H]. Arising from the Weyl- 
Lanczos equations is the hnear tensor wave equation 

^Labc + '^Rc '^abs~ I^a ^bcs~ Rb ^cas— dacR^^ Lif,s+ ghcR.'''^ Lias — —RLabc = Jabc^ 

(8) 

where ^ ^ 

Jabc = 2^c[a;b] " Q9c[aR;b] (9) 

and 

^Labc = g^"^Labc-s;m , (10) 

from which Penrose's non-hnear wave equation^ for the spacetime Weyl ten- 
sor 

R 

^Cabcd - Cab ^^Csmcd + asm[cC"^df + -^Cabcd = </[afe][c;d] + <^[cd][a;6](ll) 

~ '^[ab][c;d] ~ "^[cdllaib] 

was derived in jl5j . 

It was shown by Bampi and Cavigha that the 5 dimensional version 
of equation (P) with the appropriate equations generalising @ , © , ® , © 
constitutes a system of linear first order partial differential equations which 
is also an EDS in involution. In this paper we are going to show that like 
the Weyl-Lanczos problem the spacetime Lanczos tensor wave equation 
can be rewritten as an EDS in involution, a result which is confirmed 
using Janet-Riquier theory. 



II. The Weyl-Lanczos Equations in 4 Dimensions 
as a System in Involution 

In ^ it was shown that the Weyl-Lanczos relations are a system in involution 
with respect to the spacetime variables. The Weyl-Lanczos equations always 
constitute of a system in involution as opposed to the Riemann-Lanczos 
equations, even for vacuum spacetimes when Cabcd = Rabcd, because each 
problem is based on different equations. 

A. The Weyl-Lanczos Equations as an EDS 

The theory of exterior differential systems (EDS) can be found in many 
places such as in jH 121 1201 113j . A review together with some results on the 
Riemann-Lanczos problem in 4 dimensions is given in 0. 

The EDS for the Weyl-Lanczos equations, which form a Pfaffian system, 
was already given in p^. The Weyl-Lanczos equations together with the 
differential gauge condition © form a system of 16 linear first-order PDEs. 
In a local coordinate frame, these 16 equations look like 

fabcd Cabcd Rabcd ~l~ Rabdc Rcdab ~l~ Rcdba ~l~ ^ adi^nbc ~l~ Rncb) 

~^aciRnbd + Lndb) + ^bci^nad + L^da) ~ ^bdi^nac + Lnca) 

J-n n^-^ P J- n P — n n"""^ P — n P 

lybciJ ^nads ~r yadiJ ^nbcs ybdy ^nacs yacy ^nbds 

— n n^^(V"^ T J- r™ T J- T \ — n n^^(V"^ T 

ybcy \^ ns^mad ' ^ as^nmd ' ^ ds nam) yady ns^mbc 

bs nmc ~r i f,g-i^nbm) ' ybdy \^ ns^mac ~r -L as^nmc ~r i ^gi^nam) 

'^QacQ^^ {^^sRmbd + ^bsLnmd + ^dlLnbm) ■, 

fab — 9 {Pabns ^ as^mbn ^bsRamn ^ns^abm) ; (1^) 
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where fated = denotes the Weyl-Lanczos equations and fat = Lab^.g = 
the differential gauge conditions in local coordinates. When we construct 
the corresponding EDS, we introduce the local coordinates [x*^ , Labc, Pabcd) 
on the jet bundle J'^(M^, M^^) which form our formal manifold A4 of formal 
dimension 7V = 4 + 16 + 64 = 84. 

The exterior derivatives of all equations in (|12|) constitute our first 16 
one-forms. We also have to add 16 contact conditions K^bc in order 
to make sure that the Pabcd can be considered as if they were the partial 
derivatives of the Labc- These additional 16 one-forms are locally given by 
Kabc = dLabc — Pabcedx^ . In this Way we obtain the Pfaffian system V 

dfabcd = [Cabcd,e + O-abcde + labcde\dx'^ — dPabcd + dPabdc — dPcdab + dPcdba 
+gbc9'''dPnads + 9ad9 dPnbcs 9bd9 dPnacs 9ac9 dPnbds 
~^^adidLnbc + dLncb) — ^aci^Lnbd + dLndb) + ^bci^Lnad + dLnda) 
~^bdi'^-^nac + dLnca) ~ 9bc9^'^ {^Tis'^^mad + ^^s'^^nmd + ^^s'^^nam) 

~9ad9^^ (^ris^^nibc + ^bs^^nmc + ^^dLnbm) + 9bd9^^ (^^s^Lmac 

jSpm JT J_V'"^rlT \ J- n n'^^ (V"^ rl T J- V"^ rl T 

as^'^nmc i -L cs'-''^nam ) i yac9 ns'-''^mbd ~r i ^gUi^nmd 

'^^dl'^^nbm) J 

dfab — i,Pabns9 ,e -^mbnJ^ as,e Lamn^bs e -^ahmX ns^e)^"^ 
-\-g {dPabns ^as^^mbn ^ bs^^imn ^ns^^o.bm)i 
Kabc = dLabc ~ Pabcedx^ , (13) 

where aabcde and 'jabcde are given in Appendix A. Now, a Vessiot vector field 
for (fT^ . where 16 of the totally 64 Vabcd are determined through dfabcdiV) = 
0, dfab{V) = from (O, wih be of the form 

V = V^^^ +V^PabceJ^ + Vabcd (U) 
Ox"" OLabc OPabcd 

Because KabdV) = has to hold as well, we get Vabc = V^Pabce as used 
above. We can determine the associated system A{'P) of (|13p which is given 
locally by 

{dfabcd, dLab%, OJ^ , d Pabcd} , 

where := dx^ and where we only include 48 of the 64 dPabcd- This is 
because 16 of them can be expressed by means of the exterior derivatives 
of the Weyl-Lanczos relations and the differential gauge equations. We see 
that dim(^('P)) = 84 = c, the class of V, so that dim(P)=0, where V is 
defined as 

P = {y e V/dO'^iX^Y,) = OVX e PVl < a < s} , (15) 

where s denotes to the total number of independent 1-forms 0° in our Pfaf- 
fian system V. If we choose Y S ChariV), we have to be able to express 

Y\dKabc = Y^'dPabce - Yabcedx"" 

by means of a non-trivial linear combinations of the forms Kabc, dfab, dfabcd 
with linear multipliers A" ^ , 

^a'b'c'd' ynm 

yabc -\- ^ c d ^f^,^,^,^, _)_ \^"- dLmn ■ (16) 
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However, no such linear combination can exist because some of the dPabcd 
wiU fail to occur in one or other of the dfabcd or the dfah and in this way we 
know that no Cauchy characteristics can occur. 

Next, we want to compute the reduced Cartan characters by hand for 
this Pfaffian system V using the tableau matrix. First, we must com- 
plete {dfabcdidfab-,^^) so that it becomes a complete coframe on A^, say 
{dfabcd,dfab,Kabc,^'^,'^^), where we use the forms in (|T3|) as cobasis ele- 
ments. Accordingly, we must add the 48 new cobasis elements 

vr^ ^ dPabci, A = l,---,16, 
TT^^ dPabc2, A = 17,---,32, 
vr^^ dPabcS, A = 33, •••,48, 

where the ordering of the Pabcd based on the indices abc is given in Ap- 
pendix A and where A is a collective index subject to Einstein's summation 
convention and corresponding to the set of indices abc . Here, we write 
io^ = dx^, uP' = dx"^, Lo'^ = dx^, Lo"^ = dx^ to form the independence condi- 
tion given by 17 = W"*^ A A A uj^. Using this notation we can express the 
16 exterior derivatives of the contact conditions, where the a are arranged 
in the same way as the Labc in Appendix A, as 

de"" = Al^TX^ A + Al^TX^ A + A^^^^ Aoj^ - dPabcA A . (17) 

Further calculations are based on the assumption that we can express each of 
the 16 dPabci in (fTTj) as a distinct linear combination of the dfabcd, Kabc, dL^^b ■« 
and the oj^. Later on, it will be easy to verify that this is true for all space- 
times with diagonal metric tensor. The tableau matrices A^2)^A3 ^A4 
can now easily be determined from ()17() . where we obtain the only non- 
vanishing components to be — ^as ~ ^A4 ~ ~^ ^ ^ refer to 
the same group of indices abc. We do not need to determine AJ^^ explicitly 
as it does not contribute to the characters. This is because the sum of the 
ranks of the tableau matrices of increasing order cannot exceed 48 and so 
the terms in AJ^^ do not contribute to the rank. In this way, the matrix 

I Ki \ 

AO 

^A2 
^A3 
V ^A4 / 

is a 64x48-matrix which looks like: 

/AO \ 

0^0 

.4 
V Xi X2 Xs J 

where Xi, X2, X^ stand for A'^^ and each A is given by the 16xl6-matrix 
A := — Ii6- From this we can immediately deduce that the reduced char- 
acters are s'l = s'2 = S3 = 16 and s'^ = 0. However, this method does 
not supply Sg but here clearly Sg = s = 32 which is simply the number 
of independent 1-forms in (|13)) consisting of the 10 exterior derivatives of 
the Weyl-Lanczos equations, the 6 exterior derivatives of the differential 
gauge conditions and the 16 contact conditions Kabc- The reduced char- 
acters (sq, s'j, §2, S3, S4) then are given by (32,16,16,16,0). When (fTT?|) is 
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pulled back onto the submanifold, where fabcd = and fab = 0, we get 
(16,16,16,16,0). 

Lastly, we make an Ansatz for a transformation which absorbs the re- 
maining torsion terms -B^- . Torsion terms can arise when we write each d6" 
of a Pfaffian system V as 

= A^.TT^ A u;^ + A + ^C^y A vr'^ . (18) 

The only non-vanishing torsion terms here are those of type -B^^. In order 
to make them vanish, we must look for a transformation <I> of the form 

where we choose the in such a way that Bf^ = 0. Such p^ must be 
solutions to 

d = Bf^-Al,p\, (19) 

where we can choose p\ = Q because the Bf^ always vanish due to skew 
symmetry. For any such solution of the pf the torsion is absorbed and the 
system is in involution. This rather cumbersome calculation can be carried 
out using a REDUCE code in [Hj based on the package EDS P]. We obtained 
the result that the Cartan characters coincide with their reduced counter- 
parts and that the torsion can be absorbed for a number of spacetimes such 
as Schwarzschild, Kasner, the Debever-Hubaut class, Godel, the pp-wave 
spacetimes and conformally flat spacetimes. 



B. The Weyl-Lanczos Equations as a System of PDEs 

The theory in this section and for the Lanczos wave equation below is based 
on a modernised version of Janet-Riquier theory jlUl 1111 117j . In this mod- 
ernised form Janet-Riquier theory can be found in |16[ [T!1] and, a review 
together with some results on the Riemann-Lanczos problems in 2 and 3 
dimensions is given in |H1E]. 

The 10 independent Weyl-Lanczos relations and the 6 differential gauge 
conditions in local coordinates are given by (|12j) . Using the first computer 
code in Appendix C in 0, we can derive the symbol M.i for any spacetime 
with diagonal metric 

ds^ = ai{dx^f - a2{dx^f - a^idx^f - ai{dx'^f , (20) 

where 01,02,03,04 depend on all spacetime variables. In order to obtain a 
ranking, we replace the 4 components L121, -^^131, -^>i4i, -^^122 and their partial 
derivatives by solving @ for them and we choose an ordering >- for the 
Pabce in such a way that Pabd >- PabcS ^ Pabc2 >- Pabci and then the sets 
of Pabce Ordered according to TZ^'^^ for each e = 1,2,3,4, where TZ^'^^ is 
given in Appendix A. This produces an orderly ranking and induces such 
a ranking amongst the symbol variables Vabcd |19| I16j. where we can solve 
each equation for a different variable Vabd corresponding to a Pabci- The 
number in brackets indicates the corresponding Weyl-Lanczos equation or 
differential gauge condition given in Appendix A. Then, M.i is given in 
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orthonomic form by 
n~| ^3434 = 



1/1331 H n441 












m 



13 



14 



15 



16 



^2434 
^2424 
^2334 
^2324 
^1424 
^1434 
^1334 

Vl234 
^1324 
Vl244 

^1344 
ViM4 

V2344 
V2444 
V3444 



04.. 
1/1321 

ai 



ai 
04 



04.. 
— 1/2332 
02 



— V2442 + V3443 

^2 



V2322 + V2443 + V3442 

0-2 

04t7 , "2041 

1/2323 + ^^2442 

03 



V2343 H ^^2422 

02 



0103 
V3432 



^1431 + Vs 



02 
01 

04.. 
— 1/1323 
03 
04 



2342 



2423 



^1331 


02.. 
V3 




03 


_ 


Vu21 


01 




O2 


^3433 


03 





3443 



04 

^1233 + ^1442 

03 



04.. 

1/2331 

03 



02 



1232 



— 1/1322 + 1/1443 H 1/2321 

02 02 



V1343 H ^^1422 

02 

^1243 + ^1342 



1/2421 



02 
^1432 ■ 



1/1242 

02 



2341 



2431 



^1243 + ^1342 — ^142; 


+ ^2431 




^^2331 — 


^2441 + 
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03 


04.. 
1/1322 — 
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— V3441 
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03 
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04.. 
1/1422 — 


04.. 
1/1433 


H 1^2421 
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03 


02 


03 


1/1231 + 


04.. 
1/1321 


04.. 
1/2322 


04.. 
V^2333 


ai 


01 


02 


03 


1/1241 + 


04.. 
V^1421 


04.. 
1^2422 


04 

K2433 


01 


01 


02 


03 


1/1341 + 


1/1431 


H 1/2342 


04 

^2432 


ai 


ai 


02 


02 


^V3433 . 








03 









^1233 



(211 



All 16 equations in ()21() are of class 4 because all variables x^, 
multiplicative variables for each equation so that we obtain = 0, (3^^ 



0, P]; ' = 0, Pl"' = 16, where the upper index (k) in Pq indicates the num- 
ber of multiplicative variables and the lower index the order of the system 
of partial differential equations. For the sum X^^^^ Pi'^ = 16, and because 
we only have 16 equations, this means that we are already using coordinates 
which are 5-regular. One must use 5-regular coordinates, which are coordi- 
nates which gradually maximise Pq"'\ then Pq"" +Pq^^ and so on, to ensure 
the results obtained are intrinsic. 
Then, the Cartan characters are a 



,(4) 



(1) 



— 16, 

They are computed according to the formula 



(2) 



16, a 



(3) 



16, a 



(4) 



a 



(fc) 



m 



n + q 



k 
1 



1 



/3f 



0. 



(22) 



of which details can be found in jl9l 1161 |S1 1^. Now, we must verify that 
our symbol is involutive. Firstly, we prolong pi|) by differentiating each 
equation in TZi with respect to , x'^ , , . This leads to A42 consisting 
of 64 equations. We find that we can rewrite the prolonged symbol A42 in 
such a way that each equation contains a distinct component Vabcde given 
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in Appendix A so that r{^A2) = 64. But we also have Ylk=i ^ ' /^i^^ ~ 
which means that TWi is involutive, where further details on this calculation 
are given in [Sj. In order to show that the system is formally integrable, we 
have to verify that the canonical projection of 7^2 from second to first order 
'^i{T^2) coincides with TZi itself, which means that 7rf(7^2) = T^i- But our 
system is a linear system of PDEs and the only way an integrability condition 
can arise, is, if in any of the prolonged equations all the second-order partial 
derivatives Sahcde defined by Sahcde = Labc,de can be eliminated completely. 
This is impossible because r(A^2) = 64 and so formal integrability follows. 
We conclude that we have verified by using Janet-Riquier theory that the 
Weyl-Lanczos equations together with the differential gauge condition form 
a system in involution for spacetimes with diagonal metric tensors. 

If we wish to look at the symbol for a general spacetime, the calculations 
become more cumbersome. Therefore, we decide to look at the plane-wave 
limit which all spacetimes possess, and we shall analyse this system instead. 
A good account of the plane- wave limit of spacetimes is given in There, 
a part of a properly embedded null geodesic 7 is taken and a corresponding 
procedure applied which leads to W-y - a plane wave limit. All spacetimes 
can locally be expressed using a line element 

ds"^ = 2dx^dx^ + a{dx^f + Ih^dx^dx^ + 2h/^dx^dx'^ - C33{dx^f 

-2c3idx^dx^ - CAi{dx^f , (23) 

where a, 63, 64, 033,034, C44 are functions of all 4 coordinates. When a plane- 
wave limit is taken, the metric H'23() becomes jl4| 

ds'^ = 2dx^dx^ - Cssidx^f - 2Cudx^dx^ - C^idx^f , (24) 

where C33, C34 and C44 are arbitrary functions of x^ only. We can determine 
the symbol for the metric (|24|1 but we see that for whatever ranking we 
choose, the coordinates are not (5-regular. This means that we would have 
to prolong the system to second order and see whether we can obtain the 
desired result for the prolonged system. But in order to avoid this, we 
perform a linear coordinate transformation which we choose to be 

dx^ = aiidx^ + ai2dx'^ , 

dx^ = a2idx^ + a22dx'^ , 

dx^ = dx^ , 

dx^ = dx^ , (25) 

where 011,012,021)022 are arbitrary constants. To simplify matters, we first 
look at the special case on = 012 = 021 = -^,022 — which produces 

a new, only slightly different metric line element for (|24j) which is 

ds^ = {dx^f - {dx^f - C33idx^f - 2C3idx^dx^ - d^idi^f . (26) 

Even though this seems a minor transformation, the change from charac- 
teristic to non-characteristic coordinates is necessary in order to obtain 
intrinsic results for the Cartan characters because the new coordinates (x^) 
are 5-regular coordinates. Then, we evaluate ()12|) for the new line element 
(PH) and order the Pabcd such that Pabd >- Pabc2 >- PabcS >- Pabd instead 
whereas we leave the ordering in each set Pabce unchanged as in 7^^'^^ which 



8 



leads to an orderly ranking. We can now solve each of the symbol equations 
for one distinct variable Vabd and so obtain an orthonomic system 

^1211 = -^(^1222^ + ^1233(^44 — ^1234(^34 — ^1243(^34 + ^1244(^33) 

^1311 = -^(^1322^ + ^1333(^44 — V1334C34 — V1343C34 + V1344C33) 

^1411 = ■^(— ^1223(^34 (^33^ A — F1224A + Vi333C^^C34C44 — ^1334(^44 

— ^1343(^33^(^34^ + V1344C34 + F1422A + Vi433(C44 — C^^C34^)) 
^1221 = -r(^1212A — F1313C44 + V1314C34 + F1413C34 — Vi4i4C33 



A 

+^2323(^44 — ^2324(^34 ^ ^2423(^34 + ^2424(^33) 
1 

A' 



^1331 — -T (^1212 ((^33(^34^ — (^33^C44) + V1313 A — V2323C34^ 



+^2324(^33(^34 — V2332A + V2423C33C34 — ^2424(^33^) 
^2321 = -^(^1223^ — V1232A + V1322A + V1333C44 — V1334C34 — F1433C34 
+F1434C33) 

^2331 = ■^(^1222((^33 — (^33^(^44) — ^1233^ — V1332A + V1423C33C34 

— ^1424(^33^ + Vi323(C33C44 — 2C34^) + Vi324 (^33(^34) 
^2421 = (^1223 ((^34(^44 — (^33"'^(^34^) — Vi242 A + ^1333(733^(734(744 

— VisuCu + Vi343(733^A + V1422A — 14433(^33^(^34^ + V1434C734) 
^1231 = (^1213 A + V1312 A — F2322 A + V2343(^34 — ^2344(733 — ^2433(^34 

+V2434C33) 

^1321 = ■^(^1213A + F1312A + V2333(744 — V2334(734 — ^2433(^34 + ^2434(733) 

^1241 = ■^(^1214A + V1412A + V2343<^44 — ^2344(7*34 — V2422A 

+V2434((7'34 — (7*44)) 
^1421 = (^1213 ((^34(744 — (7^^(734^) + Vi412 A + ^2333(^33^ (^34 (7*44 

— V2334(744 + V2343(7^''^A — V2433(7^^(734^ + V2434(734) 

^1341 = ■^(— ^1212(734 A — Fi3i4 A — 1/2323(^33(^44 + V2324 (734^ + ^2423(^33 

C44 — V2424(733(^34 — V2432A) 
^1431 = ■^( — Vl212(734A + yi4i3A — V2323(^34(744 + ^2324(733(^44 — ^2342A 

+^2424(734^ — ^2424(733(734) 
V2341 = ■^(— Vl222(7'34A — V1234A — Fi323(734(744 + ^1324(733(744 

+^1423(^33(^44 — ^1424(733(734 — F1432A) 
V243I = ■^(— Vl222(734A — F1243A — Vi323(734(744 + ^1324(733(744 — V1342A 

+^1423(^33(^44 — ^1424(733(734) , (27) 

where A = C3SC44 — 0^4^. This system is composed of 16 equations all of 
class 4 which leads to /jf ^ = (3f^ = (5f^ = 0, /jf ^ = 16 and therefore to 
the Cartan characters = = of' = 16 and a^f' = 0. This result is 
intrinsic because our coordinates are now (5-regular. 
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Next, we prolong each of the 16 equations and obtain 64 equations which 
can be modified so that each of the corresponding symbol equations con- 
tains a distinct variable VabcU,Vabc22,Vabc33,VabcU or VabcSi and therefore 
r{A42) = 64. This agrees with Ylk=i ^ ' M'^'' ~ ^° that Aii is involu- 
tive. No integrability conditions can occur for the same reasons as for the 
system ()21() above and the system is in involution with Cartan characters 
(16,16,16,0). 

We obtain the same result for the plane wave limit of any spacetime as 
for the for the case of spacetimes with diagonal metric and we conjecture 
that the same system on spacetimes with arbitrary analytic metric will be 
in involution with Cartan characters (16, 16, 16, 0). 



III. The Lanczos Wave Equation in 4 Dimensions 

In this section, we look at the Lanczos wave equation at first as an EDS and 
then as a system of PDEs. We again determine the Cartan characters and 
show that it consists of a system in involution using both theories, a result 
which can be derived from the scalar wave equation directly. 



A. The Lanczos Wave Equation as an EDS 

We can also describe the Lanczos tensor wave equation in terms of an EDS 
on a jet bundle J7'^(M^,M^^) with formal dimension N = 244 for which we 
choose the local coordinates {x^, Labc, Pabcd, S abode) composed by 4 spacetime 
coordinates x^, 16 Labc, 64 Pabcd and 256 Sabcde- Here, again Sabcde are the 
variables corresponding to the second-order partial derivatives of the Labc 
when projected onto our spacetime manifold. 

First, we write the Lanczos wave equation ((HI in solved form and denote 
its components by Wabc- In a local coordinate frame Wabc is then given by 

w _ „sm r c _ r _ pn r _ pn j _ pn p 

abc y ['-'abcms ^ am,s^nbc ^ bm,s o-'n.c ^ cm^s^abn ^ as^nbcm 

- bs'- am 



I pn pfc r I pn pA: t _, pn pfc t pn p _, pn 

"r-*- as^ nm^kbc ~r -L as'- bm^nkc ~r -L as'- cm^ribk -l bs^ancm ~r 

J I pn -pk J I pn pfc J pn p , pn pA: j 

^knc ~r i fcgi ^^l^akc ~r r bs'- embank -L cs^abnm ~r i cs-*- am^kbn 

I pn pfc r I pn p/c r pn p j_ pn p/c r 

"r-*- cs-*- bm^akn ~r i cs'- nm-^abk J- jj^g-Tabcn i" -L an^kbc 

'^^'^s^bn^akc + ^^s^^n^abk] + '^^c^abs — Ra^bcs — Rb^cas 
—QacR^^Lihs + QbcR^'^Lias — -RLabc — Jabc ■ (28) 



In addition to the exterior derivatives of , we need to add two sets of con- 
tact conditions, Kabc and Kabcd-, when we write (pS)) as an EDS. Altogether, 
we obtain the Pfaffian system 

dWabc = dWabc 

Kabc — dLabc Pabcedx 
Kabcd = dPabcd — Sabcdedx^ , (29) 

where dWabc are the exterior derivatives of the components of the wave 
equation in solved form, which are locally given by 

dWabc = [^Labc + '^Rc^^abs — Ra^'bcs — Rb^cas — QacR^'^Llbs + dbcR^^ Lias 
^RLabc Jabc\,edx -|- g [dSabcms ^ am,sdLnbc ^bm,sdLanc 
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-pn JT -pn jp I -pn -pk jt , pn pfc jt 

cm,s"'^abn ^ as'^^nbcm ~r -•- as^ nm'-''^kbc ~r J- ^^i iy,^ui^nkc 

as^ cmdLnbk ~ ^ bs^Pancm + ^bs^amdLknc + ^bs^ nm^Lakc 

_i_"nn -pfc JT -pn jp _, -pn -pk jt , -pn -pk jt 

bs^ cm"'^ank ^ cs^-'-o.bnm "r -L cs-"- am^-'-'kbn ~r cs-"- bm"'^akn 

~^^cs^nmdl^abk ~ ^ms'^Rabcn + ^ms^andLkbc + ^ms^bn.d'Lakc 



+^Z,.x^rr,dLabk] + '^RrdLabs " Rn dLbcs " -R/f 

1 



ms-^ cn"'-"afefcj ~r '^^abs "'^bcs ^b "'^cas 

QacR^'dLibs + QbcR^'dLias - ^RdLabc ■ (30) 



A Vessiot vector field for is of the form 

d d d d 

V = V^-:—- + V^PabceTT^ 1- V^Sabcde-^ 1" VabcdeTT^ ; (31) 

ox" oLabc oPabcd O^abcde 

where 16 of the 160 Vabcde are determined by requiring that dWabdY) = 0. 
When we apply the 2-forms dKabc-, dKabcd to the two Vessiot vector fields 
and V"^, we can see that dx" f\dPabce{y^ ^ ^^) vanishes identically. Therefore, 
we only have to examine the second set dKabcd = —dx" A dSabcde when we 
form our integral elements of dimensions greater than one. 

When we start to compute the Cartan characters of (|29|). we obtain the 
values So = Sq = 96 again by counting the independent 1-forms in (I^U]) 
for the full and for the reduced system respectively. These 1-forms consist 
of 16 exterior derivatives of the Lanczos wave equation, 16 first-order con- 
tact conditions Kabc and 64 second-order contact conditions Kabcd totalling 
96 independent 1-forms. If we pull the system back onto the submanifold 
defined by Wabc = 0, we obtain sq = s'q = 80. 

The associated system of ()29() can again be determined and it is given 

by 

{dWabc, Kabc, Kabcd, , dSabcde} , 

where 16 of the totally 160 dSabcde can be expressed by solving dWabc for 
them so that we only add the 144 remaining dSabcde for which the dWabc 
are not solved for. We see that dim(^('P)) = 244 = c which means that 
dim(P) = so that again, no Cauchy characteristics can exist because for 
an y G Char{V) 

Y\dKabcd = dSabcde " Yabcdedx^ 

has to be such that its RHS can be expressed by means of a non-trivial linear 
combination of the form 

^abc J^^,^,^, _j_ 6 c d Xa'b'c'd' + ^^^'^ dWmns , (32) 

which is impossible. 

Again, we wish to determine the reduced Cartan characters using the 
tableau matrices A'^- and examine whether the torsion terms i??- of (|29|) can 
really be absorbed. Therefore, we must complete {dWabc, Kabc, K abed, 
to a coframe on our formally 244-dimensional jet-bundle. In (|29|) we have 
96 one- forms to which we add the 4 one- forms oj^ = dx^, uP' = dx"^, = 
(ix^, Lo^ = dx^ again composing the independence condition Q ^ so that we 
need to add a further 144 forms vr^ in order to obtain a complete coframe 
oi N = 244 one-forms. Again, A is a collective index again subject to 
Einstein's summation convention, this time replacing certain sets of indices 
abed in dSabcde- In order to obtain intrinsic values, we choose to replace all 
components dSabcde, where either d or e are = 4, first, and so on based on 
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being our first coordinate, then x^, then and last. The correspondence 

of TT^ ^ 

dSabcde is given in detail in Appendix B. Note that we did not 
replace the 16 dSabcii as this would lead to an excess of 1-forms in our 
coframe but we solve each dWabc for a distinct dSabcii- 

Because dKabc = (mod ("P)), we only need to consider dKabcd when 
computing the tableau matrix and the reduced Cartan characters. We can 
now write the 64 dKabcd as 

= Al2^^^^oJ'^+Al^-K^^uJ^+A%'K^^uJ^-dSabcll^^^+B'^jOJ'^J . (33) 

Here, we assumed that we can express all 16 dSabcii by means of a distinct 
linear combination of the dWabc: J^aba Kabcd and the lv^. From this we find 
that the only non-vanishing tableau matrix components A^- for « 7^ 1 are 
given by A^- = — 1 when the indices a and A correspond to the same set of 
indices abed as given in Appendix B and AJ^^ = otherwise. This leads to 
the tableau matrix 

/ AIa \ 
A" 

V ^li J 

which is a 160xl44-matrix of the form 



M : = 



1 ^ 










B 





c 










X2 


^3 / 



where Xi, X2, stand for AJ^^ and A is a 64x64-matrix of the form — 164) B 
a 48x48-matrix — 148 and C a 32x32-matrix —132- We see that even without 
the 16 rows for the entries X the maximal rank r(M) = 144 is obtained. 
This immediately leads to s'l = 64,^2 = 48, S3 = 32 and S4 = meaning that 
the set of reduced Cartan characters [s'q^ s'l^ s^-, s'^-, s'^'^s (96,64,48,32,0). 

As a last step, we again give an Ansatz for a transformation <!> by means 
of which the torsion terms B'^- can be absorbed. We perform a transforma- 
tion (p such that 

Q = Bt,-Al,p'\, i^l, (34) 

and we set p\:=\ because all Bf- for both i,j 7^ 1 vanish identically. This 
leads to Bf^ = for any such solution for the p^. Then, the torsion terms 
are absorbed and the system is in involution. Again, we do not carry out 
this longish calculation by hand but refer to results in confirming that 
the Lanczos wave equation consists of a system in involution with Cartan 
characters (96,64,48,32,0). For a number of spacetimes such as Kasner, 
Godel,Schwarzschild, conformally flat spacetimes, this result was shown us- 
ing the REDUCE code given in [Sj. 

B. The Lanczos Wave Equation as a System of PDEs 

The Lanczos wave equation constitutes of 16 second-order equations for the 
16 independent Lanczos components, where we again imposed (jSJ and Q. 
We can determine the symbol of the wave equation and from this obtain the 
Cartan characters. We know that the wave equation has the form Q. We 
denote the system formed by these 16 components of the wave equation by 
.wave} where the index "2" refers to the order of the system. Due to the 
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definition of the symbol, only the highest-order derivatives matter, which are 
second-order derivatives here. Therefore, the only terms contributing to the 
symbol M.2^wave are parts of the nLabc-terms, namely all the g^^dsdmLahc- 
The 16 equations for the symbol for an arbitrary spacetime then look like 

+2g^^VabclA + 2g''^Vabc2Z + 25'Va6c24 + 25^'K6c34) • (35) 

We can easily see that by ordering the Sabcde in such a way that SabcH >~ 

SabcM >~ Sabc33 ^ 'S'afec24 >- Sabc23 >- Sabc22 ^ SabclA >- Sabcl3 >- Sabcl2 >- 

Sabcii based on x'^ ^ ^ x'^ >- x^ and by choosing for each set of Sabcij the 
ordering ■j^^"-'"'^'*^ given in Appendix B, an orderly ranking is achieved. 

Because all equations of form (|35j) are of class 4, it means that = 16 
is the maximal value for plf^ . Because all p!^''^ = for k < A, the decreasing 

(k) 

sums of the P2 are automatically 16, which is the maximal value, so that 
the system is given in (5-regular coordinates. Therefore, the following values 
for p!^^ = 0, p!^^ = 0, = 0, /Sg^^ = 16 lead to the intrinsic results for the 
Cartan characters which are a'i^^ = 64, a'^^ = 48, a'^'^ = 32, = . 

In order to see whether the system is in involution, we determine wave- 
Each component of the wave equation prolonged to third order contributes 
to the symbol Ad^^wave with 

= g^^VabclU + g'^'^yabcl22 + g^^VabclSS + g'^'^yabclAA + '^g^'^VabclU 

+2g^^Vabcll3 + '^g^'^VabcllA + '^g'^^VabcUS + '^g'^'^VabcUA + '^g^'^VabclSA 

11 22 33 44 12 

= 5 VabcU2 + g Vafec222 + 5 Vafec233 + 5 Vabc2AA + 2^ Vabcl22 

+2g^^Vabcl23 + 2g^'^Vabcl2A + 2g^^Vabc223 + 25'^''Kfec224 + 2g^^Vabc23A 

11 22 33 44 12 

= 5 K6cll3 + g Vafec223 + 5 Kfec333 + 5 VabcSAA + 2g T4bcl23 

+2g^^Vabcl?,3 + 2g^'^Vabcl3A + 2g'^^Vabc233 + 2g'^'^Vabc23A + 2g^'^Vabc33A 

= g^^VabcUA + S'^^Kfec224 + 5'^'^Kfec334 + g'^^VabcAAA + ^g^'^VabcUA + ^g^"^ 

VabcUA + ^g^'^VabclAA + 2g'^^Vabc23A + 2^^'' 1^60244 + 2g^'^Vabc3AA ■, (36) 

where 320 distinct symbol variables Vabcdef can occur corresponding to the 
320 distinct third order partial derivatives which can occur. This leads to a 
64x320-symbol matrix 

( S ••• \ 
5 ••• 



\ ••• 5 / 

where each of the 16 matrices S have the form of the 4x20-matrix 



/3" 


2gi2 




2ff" 


9^' 


2^23 


2g24 




2.934 


























\ 





5" 










25I3 


25I4 











25^3 


25^4 


^33 


2534 


<?44 






















2512 





2513 


2gi4 










2.g23 


2^24 


g33 


2g34 


544 




V 








5" 








2g^^ 





25I3 


2514 





9'' 





25^3 


2^24 


^33 


2^34 ^44 





From this we can easily deduce that the rank of S is always 4 independently 
of the choice of the metric tensor so that ''(A^s^uiave) = 

4 • 16 = 64 holds. 

But this is equal to Ylt=i ^ ' (^2^^ = 4 • 16 = 64 so that M.2,wave is involutive. 
In order to show that 'JZ2,wave is formally integrable, we need to show that 
for the canonical projection '7T2{'Tl3,wave) = T^2,wave- But again, our system is 
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linear, and the only way an integrability condition can occur for an involutive 
symbol derived from a system of linear PDEs, is, if one of the equations in 
T^3,wave Can be modified in such a way that no third-order partial derivatives 
occur. Because r{^A3^wave) = 64 is maximal, this is not possible and the 
Lanczos tensor wave equation forms a system in involution with Cartan 
characters (64,48,32,0). 

We can also derive this result directly from the scalar wave equation 

= , (37) 

where ^' is our scalar component depending on x^,x'^,x^,x'^. For the scalar 
wave equation we use a formal manifold A4 with = 19 formal dimen- 
sions of which a local basis is given by the 4 x*^, 1 dependent variable ^, 
4 Pe and 10 S^j on our jet bundle J7'^(]R^,R). When projected onto the 
4-dimensional spacetime manifold, Pe are the first-order and Sef are the 
second-order partial derivatives of ^. Trivially, the symbol of a single equa- 
tion is always involutive, no matter what the ranking we chose, is, so that 
we obtain p!^^ = (i^^ = p!^^ = and /^g^^ = 1 based on the fact that a single 
equation is always of class n. We then obtain the characters 

^ 4 ^ Q,^2) ^ 2 ^ ^^3) ^ 2 , a^"^^ = . 

We further find that r(A^3) = 4 and therefore dim(A^3) = 16 whereas 
riJZ'i) = 5 and dim(7^3) = 30. We also find that dim(7^2) = 14 and because 
the symbol is linear it 

dim(7^^^^) = dim(7e3) - dim(X3) = 30 - 16 = 14 . 

From this we conclude that dim(7?-2) = dim(7^2^^) = 14 and, because the 
system is linear, no integrability conditions can occur and the scalar wave 
equation (|37)) consists of a system in involution with characters (4,3,2,0). 
Note that in the case of the Lanczos wave equation, we are dealing with 16 
such equations because we have 16 independent components Labc and find 
the correspondence between the two sets of characters is given by 

16 • (4,3,2,0) = (64,48,32,0) . 

Lastly, we wish to compare the characters of the Weyl-Lanczos equations 
with those of the Lanczos wave equation. Therefore, we need to prolongate 
the Weyl-Lanczos equations once to second order. But for an involutive 
system the characters of the symbol A^g+r of the prolonged system together 
with the are directly determined by ^1 

^q+r — ^ ( r — I ) "^9 ' 

i=k ^ ^ 

i=k ^ ^ 

where in our case n = 4, m = 16, g = l,r = 1, which leads to the re- 
sult (48,32,16,0) for [ol'^^ ^a)^^ ^a^^^ ^01^2') for the prolonged Weyl-Lanczos 
equations. This shows that, based on the for both systems, the Weyl- 
Lanczos equations are more restrictive than the Lanczos wave equation with 
characters (64,48,32,0). The general solution of the Weyl-Lanczos equations 
only depends on 16 arbitrary functions of 3 variables whereas the general 
solution for the Lanczos wave equations contains 32 arbitrary functions of 3 
variables. 
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Conclusion 



We obtained the Cartan characters {sq, si, S2, S3, s^) for the Weyl-Lanczos 
equations given as a Pfaffian system in involution which are (32,16,16,16,0) 
and (16,16,16,16,0) when pulled back onto the submanifold where the Weyl- 
Lanczos equations themselves vanish identically. These results were obtained 
assisted by REDUCE codes based on the EDS package. The Cartan charac- 
ters could be obtained for a diagonalized spacetime and for the plane wave 
limit taken when using Janet-Riquier theory and they are (16,16,16,0). 

For the Lanczos wave equation, we showed that it also consists of a 
Pfaffian system in involution and that its Cartan characters are given by 
(96,64,48,32,0) or by (80,64,48,32,0) when pulled back. We related the Car- 
tan characters and involutiveness of the Lanczos wave equation to the cor- 
responding results for the scalar wave equation. We also found that the 
Lanczos wave equation is less restrictive than the Weyl-Lanczos equations 
allowing for 32 arbitrary functions of 3 variables as opposed to only 16. 
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Appendix A: The Weyl-Lanczos Equations 

Here, we give the expressions for the Oabcde and the "fabcde occurring in the 
EDS for the Weyl-Lanczos equations given by ((T^ . The quantity Oabcde is 
determined through 

Oabcde = ^ad,ei^nbc + Lncb) + ^bc,ei^nad + Lnda) ~ ^ac,e{^nbd + Lndb) 

~^bd,e{-^nac + Lnca) ■ (^-1) 

The quantity 'jabcde is given by 

labcde — ybci^nadsQ ,e ~ ^madK^ ns9 ) ,e — ^nmd{l- as9 ) ,e — ^namK^ ds9 ) ,e\ 

+9ad[Pnbcs9^^e ~ Lmbc(X^s9^^) ,e — Lnmc(Xb^s9^^) ,e — Ln}jm(XTs9^^) ,e\ 

— n W n"-^ — T (V"^ r,'"-^\ — T (V"^ r,^^\ — T (V"^ n'"'^\ 1 

ybdY^nacsiJ ,e ^mac\^ ) ,e ^nmc\^ as9 ) ,e ^nam\^ cs9 /,ej 

~ 9ac\Pnbds9^^e ~ LmbdiX^s9^^\e ~ Lnmd(^bs9"^^) ,e ~ Lnbm(^dl9"^'^),e] 

~^9bc,e[Pnads9^'^ ~ I^niad(^^s9"^^) ~ Lnmd(^^s9^'^) ~ -^nam{^]is9"^^)] 

+9ad,e[Pnbcs9^^ ~ I^mbc{^^s9"^) ~ I-'nmciXbs9^^) ~ -^rafem(rS5'"'^)] 

ybd,e[^nacsy -^macy^ nsy ) ^nmc\'- asy ) ^nam\^ csy )\ 

-9ac,e[Pnbds9^'^ - Lrabd{^^s9^'^) " Lnmd{^b's9^'^) 

-L„,„(r2l5"^)] . (A.2) 

When introducing the collective index A for the additional coframe elements 
vr^, we ordered the Labc according to 

-^^133 < -^^144 < -^123 < -^132 < -^124 < -^^142 < -^134 < -^143 
< -^232 < -^^233 < -^^242 < -^244 < -^234 < -^^243 < -^343 < -^344 
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from which we deduce the ordering for the vr^ <-> dPatce in such a way that 
the dPabce are arranged hke the Lake above. In this way the vr^ are labelled 
according to the correspondence 



TT 



1331 



TT^'^ ^ (i-Pl332 
Vr33 ^ (iPi333 



^32 



(iP344i , 
Vr^^ ^ (iP3442 , 



When expressing the Weyl-Lanczos equations as a system of PDEs, we used 
the standard orderly ranking (if not stated otherwise explicitely) which is 
based on the ordering Pahd >- Pabc3 >~ Pabc2 >~ Pabci- Then, amongst each 
set of the Patd we order them according to 

-P344j >~ P343i >- P243i >- -^2341 >~ P2iAi >~ P2i2i >~ P233i >- P232i 
>- Pl43i >- Pl3ii >- Pu2i >- Pl2Ai >- Pl32i >- Pl23i >~ -Pl44j ^ -pL33i • 

This leads to an orderly ranking for the Pabcd for the Weyl-Lanczos equations 
denoted by T^^'"^"* . 

Sometimes, we also use the following labelling for the Weyl-Lanczos equa- 
tions and the differential gauge conditions: 







^ /l212 
^ /l314 
^ /l234 



10 



^ /l213 1 
^ /l323 7 
^ /l324 



[3] ^ /l313 
^ /1223 



12 



^ /1214 
^ /1224 



13 



L 



14 



14 



^23% , ^ ^ ^24 , ^ ^ L34 



Lastly, we give independent symbol components Vabcde for each of the 64 
equations for the prolonged symbol M2 for the Weyl-Lanczos equations 
which are: 





i = 1 


i = 2 


i = 3 


i = 4 


9j/l212 


^14411 


V23322 


V34433 


V34344 


difl213 


^13211 


V34422 


V24433 


V24344 


dif 1313 


^13311 


V24422 


V23233 


V24244 


dif 1214 


^14211 


V!34322 


V23433 


V23344 


dif 1314: 


^14311 


^23422 


V24233 


V23244 


dif 1323 


V233II 


^14422 


^13233 


^14244 


dif 1223 


V232II 


V12322 


V14433 


V14344 


dif 1224 


V242II 


^12422 


^13433 


^13344 


dif 1234 


V234II 


^14322 


^12433 


^12344 


dif 1324 


V243II 


^13422 


Vl4233 


^13244 


diL^2' 


s 


V244II 


^13322 


^12333 


^12444 


diLis' 


s 


V344II 


^13222 


^13333 


^13444 


diLii 


s 


V343II 


V14222 


V14333 


V14444 


diL23' 


s 


^12311 


^23222 


V23333 


V23444 


diL24^ 


s 


^12411 


V24222 


V24333 


V24444 


diLs4^ 


s 


^13411 


^24322 


V34333 


^34444 



Based on this we can show that r{A42) 
is given in jSj- 



64, where the detailed calculation 
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Appendix B: The Lanczos Wave Equation 

We also need to specify the correspondence between the vr'^ and the dSabcde 
for the Lanczos wave equation. For both i,j fixed we assume that the dSabdj 
are ordered according to the ordering of the Labc in Appendix A. For the 
first 64 TT^, where we have used the fact that x'^ >- )^ >- for 
our independent variables in order to obtain intrinsic results, we obtain the 
labelling (for the pairs of indices ij = 44, 34, 24, 14) 

vr^ ^ (i5'i3344 , • ■ ■ , vr^^ ^ ^534444 , 

Tt'^ ^ (i5'i3334 , • • • , vr^^ ^ ^534434 , 
vr^^ ^ ^513324 , • • • , TT^^ ^ d534424 , 
■K^'^ ^ d5'i33i4 , • • • , TT^'' ^ d5344i4 . 

Then, the next 48 vr^ are given by (for ij = 33, 23, 13) 



7r65 ^ 


'^5'l3333 , ■ 


• • , 7r^° ^ C?534433 , 




'^5'l3323 ) ■ 


• ■ , Tl^^ ^ (iS'34423 , 


7r97^ 


'^5'l3313 ; • 


• • , 7r^^2 ^ ^534413 



The last 32 vr are denoted as follows (for ij = 22, 12) 

7rll3 ^ dSi3322 , • ■ ■ , Vrl28 ^ ^^34422 , 
7r^29 ^ (iS'i33i2 , • • • , TT^^^ ^ ^£'34412 • 

We assumed that all the dSabcii can be solved for by the 16 exterior deriva- 
tives of the wave equation dWabc and expressed the dSabcii in this way. 

An orderly ranking -/^J^"^'^'^ for the 160 second-order partial derivatives 
Sabcde for the Lanczos wave equation is based on the ordering 

SabcAA >- SabcM >- Sabc33 >~ Sabc2A >~ Sabc23 >~ Sabc22 
>- SabclA >- Sabcl3 >~ SabcU ^ Sabcll 

and then amongst each set Sabdj for ^jj fixed in the same way as for the 
Pabci in Appendix A, namely: 

5344jj >- S'343jj >~ S2A3ij >- S23Aij ^ 5'244ii ^ 5'242jj >- S233ij >~ S232ij 
>- Susij >~ Siuij >- Si42ij >- Si2Aij >~ Si32ij >- Si23ij >- Simj >~ Sissij . 



The Penrose wave equation for Cabcd can be found for all dimensions 
n < 4. It was not derived from the Lanczos wave equation, which was given 
for n = 5 by Edgar and Hoglund , and we postpone the discussion of the 
5-dimensional Weyl-Lanczos problem and the related tensor wave equation 
to another paper. 
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